Adaptively measuring the temporal shape of ultrashort single photons for 
higher-dimensional quantum information processing 
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A photon is the single excitation of a particular spatiotemporal mode of the electromagnetic field. 
A precise knowledge of the mode structure is therefore essential for its processing and detection, 
as well as for applying generic quantum light states to novel technologies. Here we demonstrate 
an adaptive scheme for reconstructing the arbitrary amplitude and phase spectro-temporal profile 
of an ultrashort single-photon pulse. The method combines techniques from the fields of ultrafast 
coherent control and quantum optics to map the mode of a fragile quantum state onto that of an 
intense coherent field. In addition, we show that the possibility of generating and detecting quantum 
states in multiple spectro-temporal modes may serve as a basis for encoding qubits (and qudits) 
into single, broadband, ultrashort, photons. Providing access to a much larger Hilbert space, this 
scheme may boost the capacity of current quantum information protocols. 



Most of the experiments performed so far in the field of 
quantum optics have relied on the generation, manipula- 
tion, and detection, of single photons and other nonclas- 
sical light states in single, well-defined modes [1 . Gener- 
ation of quantum light in a superposition of two spatial or 
spectro-temporal (S/T) modes has been used to encode 
quantum information (e.g. time-bin schemes [2 , 3 ) or for 
metrological applications [4]-[7]. However, simple single- 
or two-mode systems limit the capacity in the commu- 
nication, manipulation and storage of quantum informa- 
tion and, equally importantly, are not able to seize the 
complexity of realistic states in the laboratory. 

On one hand, the possibility of using multimode quan- 
tum states of light, i.e. a single beam holding several 
independent quantum channels, may offer several advan- 
tages. As much as wavelength multiplexing has revolu- 
tionized the area of optical telecommunications, multi- 
mode states might boost the complexity of quantum net- 
works and enhance the execution of quantum information 
protocols [HHH] • In the spatial domain multimode states 
have been used in quantum imaging applications jTBHlT] 
and the orbital angular momentum of single photons has 
been explored [TSl [19] . In the spectral and temporal do- 
main they have been proposed for cluster-state quantum 
computing [20l [21] and for enhanced time metrology [22 . 
The encoding of quantum information in the S/T degrees 
of freedom of a single photon has also been proposed for 
novel quantum cryptographic schemes, such as differen- 
tial phase shift QKD [23 . 

On the other hand, the modes where the quantum 
states are prepared are often not perfectly coinciding with 
those used for their processing and detection, and this 
may significantly degrade the quality of all the possible 
applications based on such states. For example, future 
quantum networks require that light not only conveys 
information through optical links but also interacts with 
atomic species, allowing one to perform quantum process- 



ing and implement memory units [24l [25] . These tasks 
require a very specific and precise preparation of the pho- 
ton wavepacket, i.e. its S/T mode, such that it optimally 
couples to the different possible interfaces. This goal has 
already motivated a series of theoretical proposals [26l[27] 
and experiments, where the capability of arbitrarily shap- 
ing the phase [28^ and amplitude [29 profile of the mode 
of "long" (typically hundreds of nanosecond) single pho- 
tons has been shown. 

Even the seemingly simple task of detecting a sin- 
gle photon may be completely hopeless if an accurate 
knowledge of the exact spatiotemporal mode occupied 
by the photon is not available in advance. In a typi- 
cal experiment for the complete tomographic reconstruc- 
tion of some quantum light state, homodyne detection 
only works with sufficient efficiency if the spatiotempo- 
ral mode of the reference classical coherent beam (the 
so-called local oscillator, LO) is perfectly matched to 
the signal mode [30 . If little or no prior information 
on such signal mode is at hand, or if the mode itself has 
been somehow distorted during the propagation from the 
source to the detector, one may completely miss the tar- 
get in the detection stage. 

If we restrict ourselves to a single spatial mode and just 
examine the spectral and temporal degrees of freedom, 
a monochromatic single photon at frequency uj can be 
defined as 

|c.)=atH|0) (1) 

where a) is the bosonic creation operator. A single pho- 
ton in an arbitrary S/T mode can thus be written as 

= j dwM/H \u) (2) 

and, unless the shape of the wave-packet defined by the 
complex amplitude ^(uj) is accurately known, any char- 
acterization or use of such a quantum state is doomed to 



fail or perform very poorly. Measuring and controlling 
^(cj) is therefore essential. 

Starting from the above simple and prototypical case, 
we have developed a novel scheme that combines tech- 
niques from the field of ultrafast coherent control and 
quantum optics to completely reconstruct ^(cj), namely 
the amplitude and phase profile of an ultrashort single- 
photon, by mapping it to that of an intense coherent 
field. More in detail, the shape of the LO in a pulsed 
homodyne detection system is first iteratively adjusted 
to best match the single-photon shape by means of an 
adaptive genetic algorithm based on the measured homo- 
dyne data; successively, the spectral amplitude and phase 
of the single photon are characterized by measurements 
on the optimized local oscillator. Besides demonstrating 
that this approach allows us to perform the full recon- 
struction of quantum states in unknown and arbitrarily- 
shaped modes, we also show that it may serve as a basis 
for novel quantum information protocols based on the 
encoding of qubits and qudits onto the temporal and the 
spectral degrees of freedom of single broadband ultra- 
short photons. 

The experimental setup is schematically shown in 
Fig.l. We conditionally generate single photons by the 
process of spontaneous parametric down-conversion in a 
nonlinear crystal pumped by the second harmonic of a 
mode-locked 80-fs pulse train at 800 nm. Detection of 
a tightly filtered idler photon heralds the presence of a 
highly-pure single photon in the signal mode. An ultra- 
fast, time-domain, homodyne detector [3T1[32] is used to 
measure quadrature data using coherent LO pulses com- 
ing from the laser. 

In general, the measured state is not pure, but rather 
a mixture of single photon and vacuum, that can be ex- 
pressed as 

,?|1)(1| + (1-,?)|0)(0|. (3) 

Apart from other systematic and relatively constant fac- 
tors (like the limited efficiency in the detectors, optical 
losses, electronic noise, dark counts, and residual spatial 
mismatch), the S/T mode matching between the LO and 
the single photons has a direct impact on the measured 
efficiency 7^, which can thus be used as a merit parameter 
for optimizing the LO shape. 

We use a spatial light modulator (SLM) placed in the 
Fourier plane of a "zero dispersion (4/) line" ^3 [34] 
to shape the LO. Here, different pulse frequencies are 
mapped into different spatial positions, and the pixelated 
SLM is used to apply arbitrary amplitude and phase 
modulation to the different spectral components in re- 
sponse to applied voltages. The mapping of the single- 
photon S/T profile onto the LO pulse is done by means 
of an optimization procedure based on an evolutionary 
genetic algorithm [35l |36] . One first generates an initial 
population of random voltage profiles (the individuals). 
Each individual is sequentially applied to the SLM, thus 
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FIG. 1: Experimental setup. The output of the mode- 
locked femtosecond laser is frequency doubled in a lithium 
triborate (LBO) crystal and then used to pump spontaneous 
parametric down-conversion in a thin (300 /jm) /3-barium 
borate (BBC) crystal. The idler photons are detected by 
an avalanche silicon photodiode (APD, Perkin Elmer AQR- 
14) after narrow spectral filtering (F) and single-mode fiber 
(SMF) transmission. The signal field is conditionally mea- 
sured by a balanced homodyne detector (BHD) triggered by 
the idler counts. The local oscillator (LO) is obtained by 
splitting a part of the laser emission with a high-transmission 
beam-splitter (HT), and is spatially mode-matched to the 
conditional single-photon mode by a combination of lenses. 
The S/T matching of the LO is performed by means of a 
spatial light modulator (CRI, SLM-128-D-VM) placed in the 
Fourier plane of a "zero dispersion (4/) line". The spectral 
and temporal (intensity and phase) profiles of the LO are mea- 
sured with a spectrometer (Ocean Optics HR4000), with an 
interferometric auto correlator, and with a second-harmonic- 
generation-based FROG device (FROG-Scan, MesaPhoton- 
ics). Pump and signal shaping are obtained by the use of a 
stabilized Michelson interferometer, and a 10-cm-long block 
of BK7 glass, respectively. 



shaping the LO pulse according to its particular config- 
uration of voltage levels (the genes). Each shaped LO 
pulse is used for the homodyne detection of the single 
photon and the observed value of efficiency r]i is assigned 
as the fitness parameter of that particular i^^ individual. 
The voltage sequences of the best individuals are then 
cloned and used to create the individuals of the next 
generation by means of cross-overs and mutations (see 
Fig. 2a). The algorithm keeps producing new generations 
in an iterative way until convergence towards a steady 
optimum value r]opt is reached, corresponding to the best 
possible mode-matching between the LO and the photon 
shape. Finally, we use a spectrometer, an interferomet- 
ric autocorrelator, and a commercial FROG [341l device 
to fully analyze the S/T intensity and phase profiles of 
the optimal LO pulses, a faithful representation of the 
single-photon ones. 

The method is first applied to the heralded single pho- 
tons directly coming from the parametric crystal. The 
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ets. The first, simple, form of shaping is performed 
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FIG. 2: Searching for the photon shape, a) Scheme of 
the genetic algorithm used to match the LO spectro-temporal 
shape to that of the generated single photons, b) Typical evo- 
lution of the fitness parameter (the total homodyne efficiency 
77) across several generations towards convergence to the op- 
timal shape. The black curve corresponds to the different 
individuals, while the red one is the average efficiency in each 
generation. Convergence is obtained already after about 30- 
40 generations, then a slow decay is observed, probably due 
to deteriorating spatial mode-matching and decreasing laser 
power during the long (about 14 hours) measurement run. 
The sudden, isolated, drops in the efficiency curve are due to 
the appearance of "bad" random mutations. 



optimal LO shape obtained for such un-modulated sin- 
gle photons is shown in Fig. 3a. An almost flat spec- 
tral phase indicates the absence of any frequency chirp. 
In order to put the method to a more stringent test, 
we need to make use of arbitrarily-shaped single-photon 
pulses. Differently from other recent schemes ^28l [29] 
using "long" single photons produced by narrow-band 
atomic samples, where electronic modulators and detec- 
tors can keep up with the slow temporal evolution of 
the photon wavepacket, here both the modulation and 
the detection stages have to use completely different ap- 
proaches due to the ultrashort duration of the wavepack- 
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FIG. 3: Frequency-dispersed single photon. Experi- 
mental FROG traces (left panels) and reconstructed spec- 
tral intensity and phase profiles (right panels) for: a) opti- 
mal LO pulse matching the un-dispersed single photons, af- 
ter a full amplitude and phase optimization; b) optimal LO 
pulse matching a frequency-dispersed single-photon after a 
phase-only optimization; c) same as b), but with a full ampli- 
tude and phase optimization procedure. The FWHM spectral 
width is about 9.4 nm in a) and c). 

by passing the heralded single-photon pulses through a 
10-cm-long BK7 glass block, resulting in their temporal 
stretching and spectral chirp. We initially impose a sim- 
ple phase-only polynomial modulation to the spectrum 
of the LO. Here, it is the (few) coefficients of the poly- 
nomial expansion that constitute the genes on which the 
evolutionary algorithm operates. While the linear com- 
ponent of the spectral phase accommodates for possible 
temporal delays between the single photon and the LO, 
the main correction applied by the mask and visible in 
Fig. 3b is in the quadratic phase term, corresponding to 
the introduction of a linear frequency chirp in the LO 
pulse. In the next run we let both the spectral phase and 
amplitude of the shaped LO pulses vary in the adaptive 
algorithm. In this case a better efficiency is reached, and 
we obtain the shape of Fig. 3c, which still presents the 
same quadratic spectral phase, but a significantly differ- 
ent spectral intensity that, as expected, is much more 
similar to that of the un-modulated single-photon. The 
best detection efficiency obtained while measuring the 
dispersed single photons with the optimally-shaped LO 
pulses reaches rjopt ^ 60%, whereas using a transform- 
limited LO pulse reduces it to less than 50%. If, instead 
of a short glass block, one had propagated the single pho- 
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tons through the dispersive transmission hne constituted 
by a long optical fiber, the resulting quantum state would 
have been essentially lost to a detection system not tak- 
ing these S/T modulations into account. 

More complex modulations on the profile of the sin- 
gle photon are obtained by shaping the 400 nm pulses 
pumping the parametric down- conversion crystal. In- 
deed, it was demonstrated that the single photon con- 
ditionally generated in this way essentially inherits the 
spectral properties of the pump if the filtering of the her- 
ald idler mode is sufficiently tight [37H39] . In this case, 
we place a Michelson interferometer (MI) in the path of 
the pump pulses to the parametric crystal; the MI splits 
the pump pulse and recombines the two copies with some 
adjustable relative delay that causes a sinusoidal modu- 
lation of the spectral amplitude with a spectral period 
inversely proportional to this delay. We use two different 
configurations of the MI, with approximately the same 
delay (about 150 fs) and a mod(27r) temporal phase shift 

of or TT between the pump pulses. This results in 
placing the peak or the valley of the sinusoidal modu- 
lation in correspondence to the maximum of the pump 
spectral profile. In one case {ip = 0) the MI essentially 
acts as a spectral filter that narrows the pump spectrum, 
whereas in the other {(p = tt) it digs a hole in the mid- 
dle of the pump spectrum, thus producing a two-peaked 
configuration. In the cp = situation, the single pho- 
ton is also shown to acquire the expected narrow spec- 
trum with a relatively flat spectral phase; correspond- 
ingly, the single-photon pulse is temporally stretched (see 
Figs. 4a, b, and c). In the cp = tt case, the MI splits 
the pump into two time-delayed pulses that partially de- 
structively interfere in the overlap region. The double- 
peak temporal and spectral structure of the pump pulse 
is then transferred to the heralded single-photon, whose 
experimentally-reconstructed Wigner function is shown 
in Fig.5f, and is evident in the reconstructed spectral 
and temporal profiles of Fig.4e and f. In both cases, 
starting from average initial values smaller than 10% for 
the fitness of the individuals in the first generation of the 
algorithm, we observe a typically rapid increase of the ef- 
ficiency over the next few generations followed by a slower 
adjustment toward the optimum shape (see Fig. 2b). 

In the latter case, from both the spectral and tem- 
poral points of view, the single photon does not occupy 
just one peak or the other, but rather exists in a coherent 
superposition over both the two peaks. The coherent de- 
localization in time and frequency may thus allow one to 
re-define new spectro-temporal modes and ^2 based 
on these peaks, and therefore use them to encode, manip- 
ulate, and detect, qubit information with a single pho- 
ton. One can, for example, assign the |1)^^ state to the 
photon occupying the first spectral peak, and the ll)^^ 
state to the photon in the second one. In this perspec- 
tive, such a pump modulation allows the realization of 
the 2-i/2(|l)^^ |0)^^ + |0)^^ |1)^J single-photon qubit. 
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FIG. 4: Effects of a shaped pump pulse. Measured FROG 
traces and reconstructed single-photon intensity /phase pro- 
files vs. wavelength and time for (a,b,c) cp = and (d,e,f) 
(p = TV. The measured FWHM pulse duration and spectral 
width for the cp = case are 178 fs and 6.0 nm, respectively. 
When (p = TV we observe the expected tv phase jump in the 
center of the spectrum. In this case we also find an additional 
quadratic phase modulation that is likely connected to an ob- 
served slight change of the pump spatial mode. Correspond- 
ingly, the maximum efficiency achieved with the evolutionary 
algorithm is limited to about 45% in this situation. 



In order to prove that this is indeed the case, instead 
of just a statistical mixture of the photon occupying one 
peak or the other, we use the SLM to set the LO in the 
orthogonal spectral mode obtained by adding a 0lo = tt 
phase shift between the two spectral peaks (see Fig. 5b). 
This corresponds to performing a homodyne measure- 
ment of the state by projecting it onto the orthogonal 
2"'/'(|l)^, |0)^,-|0)^^ |l)^,)one. As expected, we find 
a vanishing efficiency in this case (see the reconstructed 
Wigner function of the corresponding vacuum state in 
Fig.5g) and, if the relative phase shift (j)LO is continu- 
ously varied, we clearly observe the corresponding cosi- 
nusoidal r] modulation (see Fig.5e), a convincing proof of 
the coherent nature of the superposition state. 

The above discussion can be readily extended to 
higher-dimensional spectral qudits. For example, in the 
current setup, by simply increasing the MI time delay 
between the two pump pulses. While the generated sin- 
gle photon is still delocalized between just two temporal 
modes, in the frequency domain it breaks up coherently 
in a series of equidistant spectral peaks with a comb- 
like structure. These distinct spectral modes correspond 
to the maxima of the sinusoidal modulation of the pump 
spectrum and their number is roughly limited by the ratio 
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FIG. 5: Probing a spectrally- and temporally- 
delocalized single photon. FROG-measured spectral in- 
tensity and phase profiles of the LO mode for: a) (j)LO — 0; 
b) (t)LO = vr. The corresponding temporal intensity profiles 
are shown in c) and d), respectively. Note the clear double- 
peak structure of the LO in the (I)lo — case and the nearly 
orthogonal temporal shape obtained in the (j)LO — condi- 
tion, e) Variation of the measured single-photon homodyne 
efficiency 77 as a function of the phase (j)LO between the two 
LO spectral peaks. Reconstructed Wigner functions of the de- 
tected state: a single photon in f) for = 0, and vacuum 
in g) for = tt. 

of pump to idler-filter bandwidth. Each spectral mode 
can be individually analyzed by our homodyne detec- 
tion scheme with an appropriately-shaped LO, and also 
their mutual coherence can be completely probed by ex- 
tending the experimental approach described above. In 
addition, one might use it to probe the intrinsic multi- 
Schmidt-mode character of quantum states naturally oc- 
curring in parametric down-conversion pumped by short 
pulses [40^-^ , which have also been proposed to support 
multimode quantum information processing [43-46 . Fi- 
nally, it should be noted that homodyne measurements 
with a shaped LO might not only be used in the detec- 
tion stage, but also in the conditional preparation of such 
multi-mode single-photon states. 

In conclusion, we have demonstrated the possibility of 
combining advanced techniques of ultrafast and quantum 
optics to measure quantum states of light in unknown 
S/T modes, in a way that may allow one to recover "un- 
readable" quantum information or pre-compensate the 
possible deformations encountered during propagation in 
arbitrary dispersive channels. Moreover, the possibility 
of generating and detecting ultrashort quantum states of 
light in multiple, arbitrary, S/T modes opens new excit- 
ing opportunities by providing access to a much larger 
Hilbert space for encoding, manipulating, and decoding, 
quantum information. 
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